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Abstract 

We study the chiral critical behaviors of QED by using Non-Perturbative Renormalization 
Group (NPRG). Taking account of the non-ladder contributions, our flow equations are 
free from the gauge parameter (a) dependence. We clarify the chiral phase structure, 
and calculate the anomalous dimension of ipijj, which is enhanced compared to the ladder 
approximation. We find that the cutoff scheme dependence of the physical results is very 
small. 
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1 Introduction 



The dynamical chiral symmetry breaking in gauge theories is one of the most important 
subjects in particle physics. Many efforts have been devoted to this problem particularly by 
solving the Schwinger- Dyson (SD) equation for the fermion self-energy mainly in the ladder 
approximation with the Landau gauge[0. Adding the four-fermi interactions, the chiral 
phase diagram and the critical behaviors have been investigated]^] and its characteristic 
structures with the large anomalous dimension of ipip operator has been applied to various 
models beyond the standard model [0, ^. In QCD, the ladder SD with the running gauge 
coupling constant, the improved ladder, gives good results even quantitatively 0. However 
the ladder SD is unsatisfactory in many aspects including the strong gauge dependence]^], 
the difficulty to proceed beyond the ladder]^, no firm base for the improved ladder. 

Quite recently, numerical analysis with the Non-Perturbative Renormalization Group 
(NPRG) p, ^ draws much attention as a promising new method to study non-perturbative 
phenomena of the field theory. One of the essential features of the method is that it does 
not employ any series expansion and therefore is expected to give converging results when 
improving the approximation, which shows a great contrast against other non-perturbative 



expansion methods such as or e expansions |jTO|, |TT 



Here we briefiy explain the general aspects of the NPRG approach. The RG equation is 
the nonlinear functional differential equation for the Wilsonian effective action, the exact 
form of which is explicitly written down. To approximate the RG equation, we project 
it onto a small dimensional subspace of the original full theory space. To improve the 
approximation, we enlarge the subspace step by step. As the lowest order, we take the 
local potential approximation (LPA)[jl2|. The sharp cutoff LPA NPRG has been applied 
to the dynamical chiral symmetry breaking in gauge theories and it successfully solves 



the critical behaviors [|T3|, jl^. When we further artificially truncate the LPA down to its 
'ladder parts', it reproduces the Landau gauge ladder SD results, even the improved ladder 
results when we adopt the running gauge coupling constant. Thus the LPA NPRG should 
be regarded as a non-ladder improvement of the ladder SD results, in the course of the 
systematic approximation of NPRG. 

In this article we proceed beyond the LPA analysis of the dynamical chiral symmetry 
breaking in gauge theories including anomalous dimensions of fields, which brings about 
two new issues to be carefully dealt with. First, we have to adopt a smooth cutoff scheme 
to avoid the sharp cutoff singularities for derivative interactions. What we call the cutoff 
scheme is the profile of a function to suppress the propagation of the low energy modes in 
the internal lines and is introduced through the momentum dependent mass terms in the 
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Lagrangian, 

1 f d^q 



2 7 (27r)4.^-q)-A^-c(AVg')-</'(q), (1) 

where function C defines the cutoff scheme. In the full space NPRG, physical results 
do not depend on the cutoff scheme. However the subspace approximation breaks this 
independence [^], and the cutoff scheme dependence must be checked explicitly. 

Second, the NPRG is not manifestly compatible with the gauge symmetries due to the 
momentum cutoff. Therefore we have to introduce the gauge non-invariant operators to 
our subspace and constrain them by certain identity, so-called 'Modified Ward-Takahashi- 
Slavnov- Taylor Identity' []TB[. This identity ensures the gauge invariance of the total solu- 
tions of NPRG. 

In the following sections, we formulate the NPRG flow equations beyond the LPA and 
calculate critical behaviors in QED. Further we investigate the cutoff scheme dependence 
and the gauge dependence in the covariant a-gauge. These first results beyond the LPA 
may give us a clue to the total reliability of the NPRG results. 

As mentioned above, the ladder SD results have been reproduced in the LPA by fur- 
ther truncating the flow equations. Note that this truncation cannot be regarded as any 
projection of the RG equation to a subspace. This is crucial for the gauge independence. 
Therefore the ladder truncation necessarily breaks the gauge invariance. In our system, 
e.g., the crossed ladder diagrams are included as well as the usual ladder diagrams, keeping 
the gauge independence of the flows (in some special cutoff scheme), while the ladder SD 
takes only the ladder parts, and thus suffers from disastrous gauge dependence M. 



2 RG flow equation for the massless QED 

We start with the evolution equation for the effective action Pa [A^, ■?/','?/'] 0, where A is 
the infra-red momentum cutoff for the quantum corrections. The effective action is nothing 
but the one particle irreducible parts of the Wilsonian effective action. The infrared cutoff 
is introduced by adding the following terms to the bare action: 

= / {^^'^(-^) ■ ■ ■ ^'^(^^ + ■ ■ ^(-q)^M^(q)} ' (2) 

where a parameter k labels the cutoff scheme and we take the values k = 1,2,3, ■■ ■. With 
larger k the cutoff profile becomes sharper and it reaches the 6'-function cutoff at = oo. 
Since the above curoff terms maintain the chiral symmetry, the effective action Pa also 
respects it. The evolution equation with this cutoff scheme is written as, 

A|^r. = + St Jc-'.(c-' + ;Ar,l.j I, (3) 
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where $ is a shorthand notation of the fields $ = {y4^,V',V'} ^ind C ^ is the matrix, 











(A/g) 





2fc+2 







(4) 



To approximate the RG equation we project it onto a finite dimensional subspace. 
Specifying a subspace defines an approximation. Here we take a subspace spanned by the 
seven chiral invariant operators, 



2 A2 



2 A2 



(5) 



where F^j/ is the field strength and m is the photon mass which compensates for the lack 
of the gauge invariance of Fa- 

Note that the chiral and parity invariant 4-fermi interactions without derivatives are 
spanned by the above two independent operators Gs,Gv. Within the non- derivative (lo- 
cal potential) approximation for multi-fermi operators, the 4-fermi fiow equations are not 
affected at all by other higher dimensional multi-fermi operators, which are composed of 
at least 8-fermi operators due to the chiral symmetry. In such case, the 4-fermi fiows 
completely determine the critical behaviors for the local potential approximated multi- 
fermi system, although higher multi-fermi operator fiows are necessary when calculating 



the dynamically generated fermion mass and the chiral condensates |T^. 

Substituting the effective action (||) to Eq.(^, we obtain the RG flow equations for the 
dimensionless coupling constants i.e. Z^,, Zi,a,m'^,e,Gs,Gv (see Fig.ll]). 
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Figure |l|: Diagrams incorporated in our approximation. The diagrams in the sohd box corre- 
spond to the ladder part and the diagrams in the dashed box are crucial for the gauge invariance. 
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The effective action actually contains two redundant parameters due to the wave 
function renormalization. We set and Z3 to be unity by rescaling fields A^, ^z^. Con- 
sequently, the RG equations read, 



Q 2 

d 



-4/i(A;)e^ 



(6) 



-\ {Shik, m") - h{k, m\ a)) - + 2e (G5 + AGy) hik), (7) 

2/ 



-2Gs - [3h{k, m') - h{k, m\ a)j e'Gs + QhiK m')e'Gs 
+2h{k)Gs{3Gs + 8Gv) - Ghik, m^)e\ 



d 



dt 

d_ 

di 



Gv = -2Gv - Qh{k, m')e'Gv - {3h{k, m') - h{k, m ,a)] e G 



'V 



+h{k)Gi^-3h{k,m')e\ 



a 



(8) 

(9) 

(10) 



where t is a cutoff scale parameter defined by A = Aoe *. There appear five independent 
functions of k, a and m^. These functions are defined by the following integrals: 



h{k) 
h{k,m^) 



iir"^ Jo 



k + 1 



dxx 



ik+2 



1 + x'^+i 

1 



(11) 



1 + X 



k+l 



1 



+ 



1 + 3;*^+^ + m?x^ 
1 



h{k,m^ 



k + l 



h{k,m^,a) 



^TT^ Jo 



P-1 



57r^ Jo 



l + x'^+^y \l + x^^^+m?x^ 



1 



(12) 



+2 



l + x 



k+l 



poo 




/ dx 


^2k+l 1 


Jo 





1 + x^^^ + m?x^ 
1 



(13) 



l+x'^+^y \1 + x^^^ / a + rn?x^ 



+ 



-X 



3A:+2 



I \^ ( 1 



1 + 0;'^+^/ \1 + x^^^ / a + im?x^ 



k+ 1 
k-1 



+ 2 



/5(/c, m^, a = 1). 



1 + x''^^ J \ 1 + + m^x^ 



(14) 
(15) 
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Here we note some important features of these RG equations and integrals Ix- The gauge 
parameter a-dependence of the RG equations comes solely from I^, which vanishes with 
the special cutoff scheme 'k = 1'. Therefore under this cutoff scheme, our RG equations 
are completely free from a-dependence. Of course this cancellation basically depends on 
the fact that our beta functions contain a gauge invariant set of diagrams, e.g., the crossed 
box as well as the box diagrams, and inclusion of the proper anomalous dimension of fields. 

The sharp cutoff (large k) limit makes the RG equations singular. The integrals /i, I2 
and I3 approach to finite functions of m^, while J4 and J5 blow up to infinity, except for 
the case of vanishing photon mass (m^ = 0). On the other hand, for = 0, integrals 
hi^k, rn? = 0) and /5(/c, = 0, a) vanish independently of k and a. The photon mass rri^ 
is to be determined as a function of the gauge coupling constant e by the requirement of 
the gauge invariance of the effective action Fa^o? and is the order of e^. Therefore, the 
a-dependence of our RG flow equations is the order of e'^. 

We should note the essential difference of our NPRG flows from the perturbative RG. 
Our RG equations describe the running of the 4-fermi and multi-fermi interactions as well, 
e.g., dGs/dt ~ G*!, which finally gives the effective potential of the composite operator 
ipip so that we can investigate the dynamical chiral symmetry breaking within our RG 
equations. Also the beta function of the gauge coupling constant e contains the so-called 
quadratically divergent diagram contribution, i.e., de/dt ~ eGs, which will need special 
care in relation to the gauge invariance. 



3 The critical behaviors of QED 

We investigate the Gg — sub-system described by the Eqs. (H) and (|) with the van- 
ishing photon mass approximation (m^ = 0). This system does not suffer from the gauge 
paremeter dependence, and also it has the infinite k limit which completely reproduces the 



sharp cutoff LPA (Wegner- Houghton equation) results [p!3|, |T^. In more precise treatment 
with the running photon mass, a small gauge dependence exists, except for the special 
cutoff profile of = 1. The gauge coupling constant is fixed, no running. Complete anal- 
yses with the running gauge coupling constant will be described in the forthcoming paper 
\L% . Our results here should be compared with the QED SD results with the fixed gauge 



coupling constant [Q]. Then the RG equations determining the critical behaviors are the 
following coupled differential equations (note l3{k,0) = l/(87r^)): 

^Gs = -2Gs + 2hik)Gsi3Gs + 8Gv) + \Gs-6l2ik,0)\^ (16) 

^Gv = -2Gv-XGv + Iiik)Gl-3l2ik,0)X^, (17) 
where A and I2 are defined by A = 3e^/47r^ and (47r^/3)^/2 respectively. 
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In the RG approach, the critical behaviors are governed by the structure of the fixed 
points. First we consider the case of vanishing gauge couphng constant, where the system 
is a pure fermi system and we are solving the Nambu-Jona-Lasinio model beyond the ladder 
(the large N leading) approximation. As is shown in Fig.|^, we have two fixed points, the 
Gaussian infrared fixed point at the origin and the modified Nambu-Jona-Lasinio (NJL) 
ultraviolet fixed point extended to the vector operator Gy 0. Also there appears another 
3rd fixed point at the negative Gq region, which we will ignore here as a fake due to the 
approximation. 




Figure |2[ The flow diagram on the Gs — Gv plane with the vanishing gauge couphng constant 
and with the cutoff scheme 'fc = 1'. 



There are two phases, the strong phase and the weak phase, and the renormalized trajec- 
tory is the straight line defined by G5 = SGy. In the strong phase, the 4-fermi interactions 
grow up and diverge at the finite cutoff scale in this truncation. This divergence itself 
may be regarded as a signal of the limitation of our truncation, and we expect that in 
an enlarged subspace this singularity will disappear. However, in the local potential ap- 
proximation, we can still estimate the dynamical fermion mass and the chiral condensate 
by taking account of the infinite number of higher dimensional fermi operators. In fact, 
investigating the effective potential of the composite operator ipip in either phases, we can 
conclude that in the strong phase the chiral symmetry is spontaneously broken while in 



the weak phase the chiral symmetry is respected [14 



Switching on the gauge interactions (Fig.^, these two fixed points move closer to each 
other and finally meet at some A(= A,,) to pair-annihilate. Beyond this critical gauge 
coupling constant A^,, no fixed point exists. The Gaussian fixed point at A = moves to 

•^To be precise, without gauge interactions, there is no uhraviolet fixed point since the 8-fermi operators 
do not have a fixed point solution. 
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some finite values of Gs,Gy, wliicli indicates tlie scale invariant infrared effective 4-fermi 
interactions generated by the gauge exchange interactions. 



Gv/47i^ 




Figure ||: The move of the fixed points with the various cutoff schemes. The ladder part LPA 
results |l| giving the Landau gauge ladder SD equivalents are ploted by the dotted line. 



We now understand the total phase structure. For the region < A < A^, there are two 
fixed points and two phases divided by the ultraviolet fixed point, the modified NJL fixed 
point. For A > A^, no fixed point appears and the whole theory space belongs to the strong 
phase of the dynamical chiral symmetry breaking. These structures generate the peculiar 
phase diagram on the A — Gg plane. We show the k = 1 gauge independent results in Fig.^, 
where the arrows show the RG flows. 



Gs/47t 



0.50 



0.25 



— The Phase Boundary 

Broken Phase 



0.00 




Figure ^: The phase structure on the renormalized trajectory with the cutoff scheme 'fc = 1'. 
We project it onto the X — Gs plane. Note that this is not equal to the section between the 
critical surface and the A — Gg plane which will be shown in FigJ^ 
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The critical surface obtained in various cutoff schemes are drawn in Fig.^ In Figs.^ and 
we see that the fixed point position and the critical surface depend largely on the cutoff 
scheme. This is seen in the RG equations where the criticality is directly determined by 
the cutoff scheme dependent coefficients (Ji and 12). For example, the non-trivial fixed 
point on the A = plane is given by G*g = SGy = l/(4/i(A;)). It should be mentioned that 
these quantities of the position of the criticality are not the physical quantities at all, that 
is, they can not be measured. Hence these large cutoff scheme dependences do not cause 
any difficulty in our approach. 




Figure ||: The cutoff scheme dependence of the critical surface in the (GsjCy, A) space. We 
show the cross sections of the critical surface with Gs — Gy and A — Gs planes. The dotted line 
is the result of the Landau gauge ladder SD equation. 



Our chiral phase criticality qualitatively coincides with the Landau gauge ladder SD 
results 0. The quantitative differences come out of two sources. Our results take account 
of the non-ladder contributions as well as the ladder parts, which guarantees the gauge 
independence. Also our cutoff profile is smooth while the SD uses the 6'-function sharp 
cutoff. Taking a large k, say k = 100, then our results are almost equal to the LPA sharp 



cutoff results within the line width in every figures. 



4 Critical exponent and anomalous dimension 

Here, we are going to evaluate the physical quantities, the critical exponent of the ultra- 
violet fixed point (z/4), and the anomalous dimension of the fermion mass operator (7m)- 
The critical exponent z/4 is obtained by linearizing the RG flow equations around the ul- 
traviolet fixed point, and it is related to the anomalous dimension of the chiral invariant 
4-fermi operators 74 through z/4 + 2 = 74. On the other hand, the fermion mass operator 
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m{'il)il) explicitly breaks the chiral symmetry, and the vanishing mass point is its ultraviolet 
fixed point. Thus we may evaluate its anomalous dimension 7^ from the linearized RG 
equation for the operator irttipil), 

d 

— mf = (1 + 7^) mf + 0(mf) . (18) 

The anomalous dimension 7^ is evaluated in terms of the gauge coupling constant A and 
the critical scalar 4-fermi interaction G'*(A) as 

7™(A) = ^ + 8/i(A:)G^(A), (19) 
where the corresponding diagrams are shown in Fig.^. 




Figure 0: Diagrams contributing to the anomalous dimension of the fermion mass operator ipif). 




In Fig.^ we plot the critical exponent z/4 and the anomalous dimension 7^ as a function 
of the gauge coupling constant A obtained in various cutoff schemes, with those obtained 
by the ladder SD equation. Although there are large scheme dependence in the position of 
criticalities, these exponent and anomalous dimension are almost independent of the cutoff 
scheme. This indicates that our method and approximation are fairly stable against the 
cutoff profile dependence. 
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0.9 
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1.000 
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0.000 


0.000 


0.000 
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(1.0000) 
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(1.0069) 


(1.0000) 


(1.0409) 


(1.0069) 



Table 1. The critical exponent and the anomalous dimension 7^. Our gauge independent 
results {k — 1) are compared with the Landau gauge ladder SD||^ and the sharp cutoff LPA[f3| 
results. 



Comparing with the ladder SD results, the anomalous dimension 7m is enhanced much, 
while the critical exponent z/4 is almost the same. This large corrections certainly come from 
the contribution of the non-ladder diagrams. Note that our results are gauge independent 
due to the inclusion of these non-ladder contributions. In the ladder approximation, these 
two quantities are not independent and satisfy a relation, 1/4 + 2 = 27^, which is given by 
the simple ladder relation of 74 = 27^, where only the scalar 4-fermi operators are taken 
into account. In our non-ladder extension, this relation is broken and generates the large 
deviation of 7^- The numerical values of our gauge independent results {k = 1) are listed 
in Table 1, comparing them with the Landau gauge ladder SD[Q] and the sharp cutoff LPA 
results. It should be noted that the value of the critical gauge coupling constant A^ 
depends on the approximation adopted. 

5 Discussions and Comments 

We solved the next to LPA non-perturbative renormalization group and obtained the 
critical behaviors of the dynamical chiral symmetry breaking in QED. Our RG system 
with a special cutoff scheme (A; = 1) exhibits the gauge parameter independence, which 
is assured by the inclusion of the non-ladder contributions like the crossed box diagrams. 
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The NPRG approach to the dynamical chiral symmetry breaking in gauge theories has 
great potentiahty. It turns out that the NPRG with our subspace approximation gives a 
beyond the ladder approximation restoring the gauge independence. It should be noted 
that the k = 1 cutoff scheme does not work if we go to the higher orders of the derivative 
expansion, since the beta function integrals do not converge. However, we expect that the 
gauge dependence should not be large since the NPRG approach may easily incorporate 
the 'gauge invariant' set of diagrams. 

We have also confirmed the stability of our method. The cutoff scheme dependences 
in physical quantities arc found to be negligible. Also wc get rather small deviations of 
physical results compared to the LPA. For example, the anomalous dimension 7^ at is 
1.170 in the LPA and 1.177 in our next to LPA respectively. In our level of approximation in 
this article (m^ = 0), this stability is actually equivalent to the cutoff scheme independence 
since the infinite k limit corresponds to the sharp cutoff LPA. 
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